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We design a circuit structure with linear depth to implement an n-qubit Toffoli gate. The proposed 
construction uses a quadratic-size circuit consists of elementary 2-qubit controlled-rotation gates 
around the x axis and uses no ancilla qubit. Circuit depth remains linear in quantum technologies 
with finite-distance interactions between qubits. The suggested construction is related to the long- 
standing construction by Barenco et al. (Phys. Rev. A, 52: 3457-3467, f995) [H, which uses a 
quadratic-size, quadratic-depth quantum circuit for an n-qubit Toffoli gate. 
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I. INTRODUCTION 

Practical implementation of multi-qubit quantum 
gates in quest of a scalable quantum computing system 
is essential. In particular, an n-qubit Toffoli gate plays 
a key role in established quantum algorithms. Examples 
include compiled circuits for modular multiplication and 
exponentiation in Shor's number-factoring algorithm 
H[ and quantum error correction codes (Bj]. For n = 3, 
the entangling Toffoli gate, which flips 'target' state con- 
ditioned on its two 'controls', is universal in reversible 
Boolean logic, see |||. Additionally, with an appropri- 
ate single-qubit gate, the 3-qubit Toffoli gate constructs 
a universal gate set for quantum computing Q. In the 
recent years, several protocols have been proposed to re- 
alize the 3-qubit Toffoli gate and its variants in different 
physical quantum technologies, e.g . , w ith superconduct- 
ing qubits Q , trapped ions |lfj| 03 > optical elements 
[12, EH]; and cavity quantum electrodynamics [LH ]. 

A common approach to implement a highly conditional 
gate is to apply decomposition which breakdowns the 
gat e into 'elementary' gates with at most one control 
|15Ml7l | . For an n-qubit Toffoli gate, this path results 
in quadratic-size, quadratic-depth quantum circuits with 
no ancilla 0, Corollary 7.6]. For the 3-qubit Toffoli gate, 
the simplest known decomposition requires five 2-qubit 
gates [l], Lemma 6.1], or exactly six CNOTs [l8| and 
several one-qubit gates. To avoid applying a long, at 
least quadratic-length, sequence of single- and 2-qubit 
gates, several methods have been proposed to directly 
realize multi-qubit gates with trapped ions [l9l [20| . neu- 
tral atoms [2l| . or superconducting qubits [22j. 

To streamline the realization of Toffoli gates condi- 
tioned on many qubits, which can speed-up the progress 
towards scalable quantum computation, both theoretical 
and experimental attempts are extremely important. In 
this paper, we propose a theoretical approach to decom- 
pose n-qubit Toffoli gates into 2-qubit gates in quadratic 
size, but linear depth, without using additional ancilla 
qubit. For this purpose, we change the usual computa- 
tional basis states |0) and |1) and propose a construction 



which exploits quantum rotation gates conditioned on 
one qubit. The proposed construction is related to the 
synthesis framework we suggested in (23| . 

The rest of this paper is organized as follows. The pro- 
posed circuit structure is introduced in Section[II] Circuit 
depth is analyzed in Section IIIII for quantum computing 
systems with arbitrary-length and finite-length interac- 
tion distance between qubits. We compare the proposed 
structure with prior constructions in Section ITVl Section 
IVl concludes the paper with further discussion. 



II. CIRCUIT STRUCTURE 

The choice of basis states in quantum computing is 
not unique and any two orthogonal unit vectors can be 
used in a 2-particle quantum computing system to serve 
as the computational basis states. Working with rota- 
tion gates R x (tt) around the x axis, we keep = |0), but 

change the other vector to 1 = R x (tt) |0) = [ —i ] T . 
Accordingly, R x (ir) works as a NOT gate which trans- 
forms to I and vice versa. Adding one and two condi- 
tions for R x (ir) constructs analogous versions of the con- 
ventional 2-qubit CNOT and 3-qubit Toffoli gates. Ac- 
cordingly, an n-qubit Toffoli gate is a 7r-rotation gate 
around the x axis with n — 1 conditionals. In cir- 
cuit diagrams throughout the paper, k consecutive gates 
with the same control lines are shown as a single gate 
with one control and k targets. Furthermore, C' l R x (8) 
is a ^-rotation gate around the x axis with i controls, 
aR x {0)b is a ^-rotation gate on b conditioned on a, and 
aiR x (6 l )(a 2 R x (9 2 )a 3 ) = a 2 i? :c (0 2 )(a 1 i? ;c (0i)a 3 ). 

Figure [T] shows a possible decomposition for a 3-qubit 
Toffoli gate. In this figure, if at least one of the first two 
qubits is 5, then the circuit applies either an identity / 
gate or R x (^ — |) = / gate to the target qubit. Other- 
wise, R x (^ + Tj) is applied which is a NOT gate. Lemma 
[T] provides a hierarchical structure for an n-qubit Toffoli 
gat e after applying an automatic synthesis framework in 
[23j . In the following, we first prove the lemma and then 
use it in our construction. 



* Corresponding author: msaeedi@usc.edu 



Lemma 1. [23|, Section 6] An n-qubit Toffoli gate can 
be implemented as (JXJ) . 
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p = a 
q=b 
r = ab{ 




p = a 
q=b 
r = ab ® c 



Figure 1. The 3-qubit Toffoli gate and its decomposition into 
2-qubit controlled-rotation gates. Two consecutive gates with 
controls on a are shown as a single gate with one control and 
two targets on b and c. 
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Proof. Equation ([1} is illustrated in Figure [5] with 
including conditional — 7r-rotation gates to restore control 
qubits. Consider the subcircuit A shown in the figure. 
Focusing on A, assume A input qubits are eij and A 
output qubits are bi for 1 < i < n — 1. Additionally, 
assume that ak is the first qubit (starting from k = 1) 
with value 0. After applying A, we have b\ = a±, bi = 
for 2 < i < k — 1, 6fc = l, and 6^ = a,; for fc + 1 < i < n. 

Now, consider the complete circuit in Figure [5] The 
case ai = is trivial because gates in A are disabled, the 
gate with control qubit ai and target qubit a n is deac- 
tivated, and other applied gates cancel out the effects of 
each other. Therefore, we assume oi = I. Note that be- 
fore applying A, each controlled-rotation gate with con- 
trol qubit ai for 2 < i < n — 1 applies n/2 n ~' 1 to qubit a n . 
Similarly, after applying A, each controlled-rotation gate 
with control qubit a* for 2 < i < n — 1 applies — 7r/2™~ 4 
to qubit a n . 

If ak (starting from k = 1) is the first qubit with 
value 0, then conditional rotation gates with controls 
a\, a.2, • • • , ak-i are activated and a #i-rotation gate with 
#1 = + + 2^ + ■■■+ is applied to the 

target qubit. However, after applying A a # 2 - r °tation 
gate with 82 = is applied which removes the effect 
of 61 given 6\ = —62- Additionally, each gate with con- 
trol qubit ai for k + l<i<n — 1 after A removes the 
effect of its corresponding gate before A. 

Finally, if a, = 1 for all 1 < i < n — 1, then all gates 
before A are enabled and all gates after A are disabled 
and a ^-rotation gate with 9 — + ^7 + 5^=^ + ■ 



7r is applied to the target qubit a n . 



□ 
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Theorem 2. An n-qubit Toffoli gate with controls 
ai,a2,--- ,a n -i and target a n can be implemented by 
a network of the form given in Figure [3] where all gates 
are conditional ^-rotation gates around the x axis. 

Proof. To prove, we restructure the circuit shown in 
Figure [3] as illustrated in Figure [5] To verify, note that 



gates in the first (top) n — 1 lines construct an (n — 1)- 
qubit Toffoli gate, gates in the first n — 2 lines construct 
an (n— 2)-qubit Toffoli gate, ■ • • , gates in the first 3 qubit 
constructs a 3-qubit Toffoli, and finally gate in the first 
2 qubits is a CNOT, these gates are specified in Figure 
|3] too. Based on Lemma [TJ the circuit shown in Figure [3] 
implements an n-qubit Toffoli gate. □ 
Figure 4(a) shows the proposed construction for a 5- 
qubit Toffoli gate. In Figure |4(b)[ the construction is 
restructured differently to better reflect the hierarchial 
structure in (23| . Comparing the decomposition in Fig- 
ure |4(by| with the conventional method in see Figure 
[SJ reveals the main differences between two methods. To 
count the number of 2-qubit gates in the proposed con- 



struction, note that there are 2S* 



"i + n — 1 gates 



to construct the transformation on the target line, and 
2S 



-n — 3 



( = 1 



i + n — 2 gates to restore control lines to their 
original values. Therefore, the total number of 2-qubit 
gates in the proposed construction is 2n 2 — 6n + 5 or 
2n 2 + 0(n). 



III. DEPTH ANALYSIS 

In this section, we show that in spite of the quadratic 
size of the proposed structure for an n-qubit Toffoli gate 
(no ancilla), circuit depth is linear. In order to consider 
depth, we restructure the construction shown in Figure[3] 
In particular, we change the structure to have gates with 
common targets (vs. common controls in Figure [3]) in 
sequence. Additionally, we divide the circuit in Figure [3] 
into 6 parts, namely C\ , C%, ■ ■ ■ , C§ as shown in the figure. 
To evaluate circuit depth, we focus on C\ . The result can 
be extended to the whole circuit. Figure [5] illustrates C\ 
in Figure 4(a) with time steps for each gate. 



Theorem 3. The proposed structure for an n-qubit Tof- 
foli gate can be implemented by a linear-depth circuit. 

Proof. Restructuring the circuit structure in Figure [3] 
to have gates with common targets in sequence, one can 
verify that in C1+C2 there are n — 1 gates with targets 
on qubit n, n — 2 gates with targets on qubit n — 1, 
• ■ ■ , one gate with target on qubit 2. Assign time steps 
1, 2, • • • , n — 1 to n — 1 gates with targets on qubit n. 
Next, consider the n — 2 gates with targets on qubit n — 1. 
Among these gates, n— 3 gates can be executed in parallel 
with the gates with targets on qubit n. Precisely, gates 
with targets on qubit n—2 can be executed in time steps 
3, 4, • ■ ■ , n — 1, n. Similarly, the next n — 4 gates can be 
executed in time steps 5,6,-- - , n+1. Following this path 
results in 2n— 3 time steps for C1+C2. Likewise, C3 can be 
parallelized to depth 2n — 5, C4+C5 can be parallelized to 
depth 2n — 5, and finally Cq can be parallelized to depth 
2n — 7. Altogether, circuit depth for an n-qubit Toffoli 
gate in the proposed construction is 8n — 20. □ 

While circuit depth in the proposed construction is lin- 
ear, our construction includes many long-distance 2-qubit 
gates. In general, restricting interactions to only linear 
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Figure 2. An n-qubit Toffoli gate can be implemented using 2-, 3- 
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(n — l)-qubit Toffoli gates based on {T}. 
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Figure 3. Circuit structure for an n-qubit Toffoli gate. Excluding gates that should be applied to restore control qubits, 
different boxes represent (n — l)-qubit Toffoli, (n — 2)-qubit Toffoli, • • • , 3-qubit Toffoli, and 2-qubit CNOT gates. These gates 
are specified by C l , for C 1 R x (tv), in the figure. The proposed construction is divided into six parts Ci,C2, ■ ■ ■ , Cs- 



dimension (ID) results in 0(n) overhead. However, cir- 
cuit depth in the proposed construction remains linear 
even in very restrictive quantum architectures with pos- 
sible interactions in a line. Assume a SWAP gate between 
qubits ai and a 2 is represented by 5(a 1 ,a 2 ). We use the 
term 'local' for gates that use neighbor qubits in a given 
architecture. 

Theorem 4. Circuit depth for an n-qubit Toffoli in 
the proposed construction is linear in architectures with 
finite-distance interactions between qubits. 

Proof. To prove, we consider ID architectures. 
One can execute a fD quantum circuit on architectures 
with interactions in a higher dimension. Working with 
C\ + C2, consider a chain of n — 1 serial SWAP gates 
S{a n , a„_i), S[a n -i,a n -2), 6>(a„_2, a.n-3), • • ■ j <5>(o2, a i) 
in sequence. For an initial qubit ordering 1, 2, • • • , n, the 
resulting ordering is n, 1,2, ••• ,n — 1 (i.e., a 1-bit ro- 
tation). Immediately after each SWAP gate, one can 
apply a local controlled-rotation gate with target on 
qubit n. Now, apply a chain of n — 2 SWAP gates 
S(a n , a n -i), S(a n -i, a.n-2), S(a>n-2,a>n-3), • • ■ j S(a^, 02) 
in sequence. Among these n — 2 gates, n — 3 gates can 
be executed in parallel with the previous gates. After 



the second SWAP chain, the resulting qubit ordering is 
n, n — 1, 1, 2, • • • , n — 2, i.e., a 2-bit rotation. Accordingly, 
we can apply n — 2 local controlled-rotation gates with 
targets on n — 1. Following this path results in 2n — 3 
time steps for SWAP gates, and 2n — 3 time steps for 
controlled-rotation gates, An — 6 2-qubit time steps in to- 
tal. Circuit size is increased by 2n — 3 for SWAPs. The 
final qubit ordering is n, n — 1, n — 2, • • • ,2,1. 

To construct a local circuit for C3 starting from qubit 
ordering n, n — 1, n — 2, • • ■ , 2, 1, we can apply the same 
structure discussed. It leads to depth 4n — 10 for C3. 
The resulting qubit ordering is 2, 3, • ■ ■ , n — 1, n, 1. At 
this time, applying the next C4 + C5 circuit is tricky be- 
cause qubit ordering has been changed from the initial 
one 1, 2, • • • , n— 1, n. Actually, the first qubit is far from 
other qubits 2, 3, • • • . For this case, we apply a linear- 
depth circuit with depth n + 5, and size 4n — 6 [H, The- 
orem 4.1] to restore the ordering 1, 2, ■ ■ • ,n — 1, n. Ac- 
cordingly, C4 + C5 , and Cq can be implemented in depth 
An — 10 and An — 14, respectively. We recover the final 
qubit ordering to the initial ordering 1,2, ••• ,n — l,n 
with another linear-depth circuit. 

Altogether, circuit depth for an n-qubit Toffoli gate 
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Figure 4. Circuit structure for a 5-qubit Toffoli gate. Circuit 
in (a) is the proposed structure. This circuit is restructured 
in (b) to use circuits for 4, 3, and 2-qubit Toffoli and CNOT 
gates. Note that direct decomposition of the gates in (b) 
results in many redundant gates. 
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Figure 5. Circuit for an n-qubit Toffoli gate in [l], Lemma 7.5] 
where V 2 = U and [/ is a NOT gate. The last gate can be 
decomposed by recursively applying the decomposition. 



with only ID interactions can be calculated as 18n — 31. 
Circuit size remains 2n 2 + 0(n). □ 
In summary, circuit depth in the proposed structure 
is only increased by a constant factor, e.g., 2.25 in ID 
architectures. Figure [7] illustrates the circuit in Figure [6] 
with only local gates. 



IV. COMPARISON WITH PRIOR ART 

The current widely-used decomposition [l], Corollary 
7.6] for an n-qubit Toffoli gate uses a quadratic-size con- 
struction with staircase structure where target of gate 
i depends on a control of gate i — 1. This results in 
a quadratic depth. The decomposition is illustrated in 
Figure In this figure, U is a NOT gate which results 
in V and where V 2 — U. The resulting multiple- 
control Toffoli gates have linear cost 48n + 0(1) in [l[ 
due to the availability of one ancilla qubit. The last gate 
can be decomposed by recursively applying the decom- 
position shown in Figure [5] using U = VNOT. Follow- 
ing this path results in controlled-ith-root-of-NOT gates 
for i = 2 1 ,2 2 ,--- , 2 n_1 . Circuit size and depth are 
48n 2 + 0(n) 2-qubit gates. 
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Figure 6. A part of the circuit shown in Figure 
tured to show parallel circuits. Numbers are t 
that gates can be executed. 
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Figure 7. Circuit in Figure [6] with only local gates. Numbers 
are time slots that gates can be executed. 



The optimizations in j25j improve the linear-cost im- 
plementation of a multiple-control Toffoli gates with one 
ancilla from A8n + O(l) to 24n + O(l). Circuit depth re- 
mains quadratic, precisely 24n 2 + 0(n). The method in 
[23l . Section 6] benefits from a recursive construction with 
quadratic-depth 2n 2 + 0(n). As discussed in Section [TT1 
and Section IlIIl our circuit size and depth are quadratic 
and linear, respectively. All methods uses gates with sim- 
ilar complexity levels for physical realization. 

In the proposed structure we assumed no ancilla qubit 
is available to facilitate circuit construction. If at least 
one ancilla exists, prior circuit structures in 0, Lemma 
7.2] and [l|, Lemma 7.3], and the extended versions (25| . 
use linear-size circuits. When 1 and n — 3 ancillae are 
available, we can apply the same circuit structures in P, 
Lemma 7.2] and p], Lemma 7.3]. Precisely, after applying 
various optimizations in (25| . we can construct circuits 
with sizes 24n — 88, and 12n — 34 if one and n — 3 ancillae 
are available — note that Peres gate has cost 4 in the pro- 
posed construction as in (25| . Reusing optimizations in 
[25| in the proposed circuit structure is straightforward. 



V. CONCLUSION 

We proposed a linear-depth quadratic-size quantum 
circuit with controlled-rotation gates around the x axis 
with no ancilla qubit to implement an n-qubit Toffoli 
gate. Restricting qubit interactions in finite length af- 
fects circuit depth and size by a constant factor. 

The proposed structure may or may not be a physically 
realizable construction in a particular quantum comput- 
ing technology. The physical implementations of quan- 
tum gates are imperfect due to various reasons including 
decoherence and error in experimental setups. In the pro- 



5 



posed circuit structure, we used ^-rotation gates around 
the x axis for 6 = ^ and 1 < k < n — 2. Obviously, 
2^=2 can be very small for large n values, which makes 
its physical implementation complicated. Small rotation 
angles may be ignored in specific applications, as done 
for approximate quantum Fourier transform |26fl . In par- 
ticular, restricting k < [log 2 n] results in e w — error. 



Since conditional Toffoli gates are building blocks for 
various quantum algorithms, in-depth characterization 
and understanding of their operations and imperfections 
possibly based on quantum tomography (2?J can be very 
useful. Recently, a multi-qubit phase gate with one con- 
trol qubit simultaneously controlling n target qubits was 
implemented using superconducting qubits [28j. Since 
we extensively benefit from such gates in the proposed 
construction, applying the method in (28| to physically 



realize conditional Toffoli gates based on the method pre- 
sented in this paper, e.g. the small circuit for a 4-qubit 
Toffoli gate, may be useful. 

Finally, while we use and I for computational basis 
states, we can also use |0) and |1). To achieve this, one 
can transform |0) , |1) to 0, f by applying n single-qubit 
gates with the same matrix M to all qubits. This should 
be followed by the proposed construction. Final quantum 
state can be restored from 0, 1 to |0) , 1 1 } by applying M'. 
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As a side note, restricting to have only one type of 
2-qubit gate can increase circuit depth and size by a con- 
stant factor given each 2-qubit gate can be implemented 
by a constant-size circuit [l|. 



A. Barenco, C. H. Bennett, R. Cleve, D. P. DiVincenzo, 
N. Margolus, P. Shor, T. Sleator, J. A. Sm olin, and 

H. Weinfurter, |Phys. Rev. A 52, 3457 (1995)] 

R. Van Meter and K. M. Itoh, 
|Phys. Rev. A 71, 052320 (2005) | 

I. L. Markov and M. Saeedi, Quant. Inf. Comput. 12, 361 
(2012). larXiv:1202. 66141 
L L. Markov and 
|Phys. Rev. A 87, 012310 (2013) 



M. 



Saeedi, 



M. A. Nielsen and I. L. Chuang, Quantum Computation 
and Quantum Information (Cambridge University Press, 
2000). 

M. Saeedi and I. L. Markov, ACM Computing Surveys 
45 (2013k larXiv:1110.2574l 

Y. Shi, Quant. Info . Comput. 3, 84 (2003), 
|arXiv:quant-ph/0205115l 

A. Fedorov, L. Steffen, M. Baur, M . P. da Silva, an d 

A. Wallraff, Nature 481, 170 (2012k larXlv:1108.3966l 
V. M. Stojanovic, A. Fedorov, A. Wallraff, and 
C. Bruder, |Phys. Rev. B 85, 054504 (2012) | 

T. Monz, K. Kim, W. Hansel, M. Riebe, A. S. Villar, 
P. Schindler, M. Chwalla, M. Henn rich, and R. Blatt, 
|Phys. Rev. Lett. 102, 040501 (200971 
M. Borrelli, L. Mazzola, M. Paternostro 
calco, |Phys. Rev. A 84, 012314 (2011)| 
T. C. Ralph, K. J. Resch, and A. Gilchrist, 
|Phys. Rev. A 75, 022313 (2007) | 

B. P. Lanyon, M. Barbieri, M. P. Almeida, T. Jennewein, 
T. C. Ralph, K. J. Resch, G. J. Pryde, J. L. O'Brien, 
A. Gilchrist, and A. G. White, Nature Physics 5, 134 
(2008). 

X.-Q. Shao, A.-D. Zhu, S. Zhang, J.-S. C hung, and K.-H. 
Yeon, |Phys. Rev. A 75, 034307 (2007)| 
J. J. Vartiainen, M. Mottonen, and M. M. Salomaa, 
|Phys. Re v Lett. 92, 177 902 (2004) | 
Shende, ST 



and S. Manis- 



V. V. 

Markov, 



Bullock, 



and I. 



IEEE Trans, on CAD 25, 1000 (2006) 



|arXiv:quant^ph/ 0406176 
M. Saeedi, M. Arabzadeh, M. Saheb Zamani, and 
M. Sedighi, Quant. Inf. Comput. 11, 0262 (2011), 



larXiv:1011.2159l 
[18] V. V. Shende and I. L. Markov, Quant. Inf. Comput. 9, 

461 (2009k larXTv:0803.2316l 
[19] X. Wang, A. S0rensen, and K. M0lmer, 

|Phys. Rev. Lett. 86, 3907 (2001)] 
[20] S! Ivanov and N. V. Vitanov, 



|Phys. Rev. A 84, 022319 (2011) 



[21] L.-M. Duan, B. Wang, 
|Phys. Rev. A 72, 032333 ( 2005) 



and H. 



Kimble, 



[22] C.-P. Yang and S. Han, |Phys. Rev. A 72, 032311 (2005)] 
[23] A. Abdollahi, M. Saeedi, and M. Pedram, Quant. Info' 

Comput. (2013k larXiv:1302. 53821 
[24] S. Kutin, D. Moulto n, and L. Smithline, Ch icago J. of 

Theor. Comput. Sci. |arXiv:quant-ph/0701194 
[25] D. Maslov, G. W. Dueck, D. M. Miller, and C. Ne- 
grevergne, IEEE Trans, on CAD of Integrated Circuits 
and Systems 27, 436 (2008), |arXiv:quant-ph/0604001| 
[26] A. Barenco, A. Ekert, K.-A. Suominen, and P. Torma, 

|Phys. Rev. A 54, 139 (1996) | 
[27] M. Riebe, K. Kim, P. Schindler, T. Monz, P. O. Schmidt, 
T. K. Korber, W. Hansel, H. Haffner, C. F. Roos, and 



[28] 



R. Blatt, |Phys. Rev. Lett. 97, 220407 (2006) | 
C.- 



Yang, 



Y.-x. Liu, 



and 



Nori. 



|Phys. Rev. A 81, 062323 (2010) 



ACKNOWLEDGEMENTS 

Authors were supported by the Intelligence Advanced 
Research Projects Activity (IARPA) via Department 
of Interior National Business Center contract number 
D11PC20165. The U.S. Government is authorized to 
reproduce and distribute reprints for Governmental pur- 
poses notwithstanding any copyright annotation thereon. 
The views and conclusions contained herein are those of 
the authors and should not be interpreted as necessarily 
representing the official policies or endorsements, either 
expressed or implied, of IARPA, DoI/NBC, or the U.S. 
Government. 



